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ABSTRACT 

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and 
the cosmological parameters which govern them. At low redshifts, the non-linear gravitational evolution 
of large scale structure produces a non-Gaussian covariance in the shear power spectrum measurements 
that affects their translation into cosmological parameters. Using the dark matter halo approach, we 
study the covariance of binned band power spectrum estimates and the four point function of the dark 
matter density field that underlies it. We compare this semi-analytic estimate to results from N-body 
numerical simulations and find good agreement. We find that for a survey out to z ~ 1, the power 
spectrum covariance increases the errors on cosmological parameters determined under the Gaussian 
assumption by about 15%. 

Subject headings: cosmology: theory — large scale structure of universe — gravitational lensing 



1. INTRODUCTION 

Weak gravitational lensing by large scale structure (LSS) 
shears the images of faint galaxies at the p ercent level and 
correl ates their measur ed elli pticitics (e.g., Blandford et a\ 
1991; pVliralda-Escudq 1 1991; [Kaiser] 1992). Though chal- 
lenging to measure, the two point correlations, and the 
power spectrum that underlies them, provide important 
cosmological information that is complementary to that 
supplied by the cosm ic microwave b ackgr ound and poten- 
tially as prec ise (e. g., Jain fc Scljak 1997 ; |Bernardeau et a| 
1997; iKaiserj 1998; |Hu fc TcgiT 3iudTl 999; IHuj 1999; |Cooray1 
199 9; Van Wacrbekc et al| 1999;~see Bartclmann fc Schnci-| 
dei 2000 for a recent review). Indeed several recent studies 
,ve provided the firot clear evidence for weak lonoing in 



ha; 

so-cal led blank fields where the large scale structure signal 



is expec ted to do minate (e.g., [ Van Wacrbekc et a j 2000 



Bacor 



et al| 2000; | Wittman et al| 2000; [Kaiser et a] 200 6} 



Given that weak gravitational lensing probes the pro- 
jected mass distribution, its statistical properties reflect 
i of the dark matter. Non lincaritico in the mono dio 



those 

tribut ion; duo to gravitational ovolution at lc 



3dahiftc 



cause s the shear field to become non Gaussian , It is well 

known that lensing induces a measureabl e three-point cor- 
relati on in the derive d convergence field ( Bernardeau et a\ 
1997; Cooray fe Hu| 2000). The same processes also in- 
duce a four point correlation. The four point correlations 
are of particular interest in that they quantify the sample 
variance and covariance of two point correlation or power 
spectrum measurements. Previous studies of the ability 
of power spectrum measurements to constrain cosmology 
have been based on a Gaussian approximation to the sam- 
ple varianc e and the assum ption that covariance is negli- 
gible (e.g., Hu & Tegmark 1999); it is of interest to test 
to what extent their inferences remain valid in the pres- 
ence of realistic non-Gaussianity. More importantly, when 
interpreting the power spectrum recovered from the next 
generation of s urveys an ac curate propagation of errors 
will be critical QHu fe White] 2000). 

Here, we present a semi-analytical estimate of the Fourier 
analogue of the four point function, i.e. the trispectrum, 
and calculate in detail the configurations that contribute 



to power spectrum covariance. Since weak lensing shear 
and convergence can be written as a simple projection of 
the dark matter density field, the problem reduces to a 
study of the trispectrum of the density field. Previous 
studies of the dark matter trispectrum have employed a 
mix of perturbation theory and non-linear scalings (e.g. 



Scoccimarro, Zaldarriaga fc Hui 1999) or N-body simula- 
tions ( [Meiksin fc White 1999). The former are not ap- 



plicable to the full range of scales and configurations of 
interest; the latter are limited by computational expense 
to a handful of realizations of cosmological models with 
modest dynamical range. 

Here, we use th e dark -matte r halo app roach t o model 



the density field ( |5eljak| 2000; flVla fc Fryj 2000a; |Scocci 



marro et al. 2000) and extend our previous tre atments of 



the two- point and thre e-point lensing statistics (Cooray et 



3- pc 

a| 2000; jCooray fc Hu| 2000). The critical ingredients are:' 



a mass function for the halo distribution, such as the Press- 



Schechter (PS; Press & Schechter 1974) mass fu nction; a 



profi le for the dark matter halo, e.g., the profile of Navarroj 
et al| (1996; NFW), and a description of halo biasing ( |Mc[ 



et al 



1997). In the mildly non-linear regime, where most 
of the contribution to lensing is expected, the accuracy of 
the halo model has been extensively tested ag ainst si mula- 
tions at th e two p oint and three point levels (Scljak 2000; 
Ma fc Fry| 2000a; [Scoccimarro et al.| 2000). We present 



tests here of the four point configurations involved in the 
power spectrum covariance. These techniques can also be 
extended to the covariance of the power spectrum of galaxy 
redshift s urveys with a prescription for assigning galaxies 



to halos dSeljak] 2000; [Bcoccimarro et al] 2000). The effect 



of non-Gaussianitics on the measured galaxy power spec- 
trum, through a measu rement of the angular cor relation 
function, is discussed in Eiscnstein & Zaldarriaga (1999). 

In §pL we study the trispectrum of the dark matter den- 
sity field unde r the halo model and test it against sim- 
ulations from (Meiksin & White 1999). In §|l we apply 
these techniques to the the weak len sing covarian ce and 
test them against the simulations of (White & Hu 1999). 



We also discuss the effect of power spectrum covariance on 
cosmological parameter estimation. 
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Fig. 1. — The dark matter power spectrum (a) and square- 
configuration trispectrum (b) broken into individual contributions 



under the halo description. The lin 



a labeled 'PTT sh 



3ws the dark 

matter power spectrum under the Peacock & Dodde (1996) non- 
linear fitting function while the curve labeled TI" is trie linear dark 
matter power spectrum (at redshift of 0). In (a), we show the power 
spectrum at redshifts of and 1. In (b), we show the square config- 
uration trispectrum (see text). In both (a) and (b), at small scales 
the single halo term dominates while at large scales halo correlations 
contribute. 



2. DARK MATTER POWER SPECTRUM COVARIANCE 

We begin by defining the power spe ctrum, trispectrum 
and power spectrum covariance in §|2.1| We then derive the 
halo model for these quantities in 



In § 2.3, we present 
results and comparisons with TV-body simulations. 

2.1. General Definitions 

The two and four point correlations of the density field 
are defined in the usual way 



(<KkiM(k 2 )) = (2^) 3 < 5 D (k 12 )P(fc 1 ). 



(1) 



(J(k0 . . . 6(k 4 )) c = (2^) 3 ( 5 D (k 1234 )T(k 1 , k 2 ,k 3 ,k 4 ) ,(2) 

where kj.j = kj + . . . + kj and <5d is the delta function 
not to be confused with the density perturbation. Note 
that the subscript c denotes the connected piece, i.e. the 
trispectrum is defined to be identically zero for a Gaussian 



field. Here and throughout, we occasionally suppress the 
redshift dependence where no confusion will arise. 

Because of the closure condition expressed by the delta 
function, the trispectrum may be viewed as a four-sided 
figure with sides k^. It can alternately be described by 
the length of the four sides hi plus the diagonals. We 
occasionally refer to elements of the trispectrum that differ 
by the length of the diagonals as different configurations 
of the trispe ctrum. 



Following Scoccimarro, Zaldarriaga & Hui (1999), we 
can relate the trispectrum to the variance of the estimator 
of the binned power spectrum 



Pi 



1 

V 



d 3 k 

v^7 



**(-k)a(k), 



(3) 



where the integral is over a shell in fc-space centered around 
hi, V S i « 4:Trkf5k is the volume of the shell and V is the 
volume of the survey. Recalling that 5(0) — > V/(2tt) 3 for 
a finite volume, 



PPj > " 



Pi) (P.. 



1 

V 



2PfS l . 



where 



Ti 



d kj f d k'j 



- T ( k.j , k^ , kj , k j 



(4) 



(5) 



Notice that though both terms scale in the same way with 
the volume of the survey, only the Gaussian piece nec- 
essarily decreases with the volume of the shell. For the 
Gaussian piece, the sampling error reduces to a simple 
root-N mode counting of independent modes in a shell. 
The trispectrum quantifies the non-independence of the 
modes both within a shell and between shells. Calculating 
the covariance matrix of the power spectrum estimates re- 
duces to averaging the elements of the trispectrum across 
configurations in the shell. It is to the subject of modeling 
the trispectrum that we now turn. 

2.2. Halo Model 

We model the power spectrum and trispectrum of the 
dark matter field under the halo approach. Here we 
present in detail the extensions re quired to mod el the 
trispectrum. We refer the reader to Cooray & Hu (2000) 
for a more in depth treatment of the ingredients. 

The halo approach models the fully non-linear dark mat- 
ter density field as a set of correlated discrete objects 
("halos") with profiles pi that for definiteness depend on 
their mass M and concentration c as in the NFW profile 
( |Navarro etaj| 1996)0 



p(x) = ^ /°h( x ~ M i> °i) > 



(6) 



and so a density fluctuation in Fourier space 

<5(k) = ^ e ' ik - x '<5 h (k;M; C ) (7) 



n 1 e l ^d h (k,M 1 ,c 1 ) 

Vi.A/i.ci 



(8) 



x This prescription can be generalized for more complicated halo 
profiles in the obvious way. 
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where we have divided space up into volumes 8V suffi- 
ciently small such that th ey conta in only one halo ni = 
n\ = rii — 1 or following Peebles (1980). The final ingre- 
dient is that the halos themselves are taken to be biased 
tracers of the linear density field (denoted PT) such that 
their number density fluctuates as 



d n 
dMdc 



(x) 



d n 
dMdc 



1 



[&o + 6i(M)V(x) + -& 2 (M)<&(x) 



( 9 ) 

where bo = 1 and the halo bias parameters are given in 



Mo et aL (1997). Thus 



(nin 2 ) 



(nin 2 n 3 ) 



-SMiSci 



dMdc 
(m) 6 12 + (ni) (n 2 ) [6q 
x (5 pt (xi)5pt(x 2 ))] ■ 



6i(Mi)6i(M 2 



(10) 



(11) 



The derivation of the higher point functions in Fourier 
space is now a straightforward but tedious exercise in al- 
gebra. The Fourier transforms inherent in eqn. (^) con- 
vert the correlation functions in eqn. ( |Tl] ) into the power 
spectrum, bispectrum, trispectrum, etc., of perturbation 
theory. 

Replacing sums with integrals, we obtain expressions 
based on the general integral 



dM / dc 



d n 
'dMdc 



bp{M) 



8 h (k u M,c)...5 h (k^M lC ). (12) 



The index /i represents the number of points taken to be 



in the same halo such that (rv{ 

The power spectrum under the halo model becomes ( Sel- 
00) 



J a 



P(k) = P lh (k)+P 2h (k), 
P lh (k)=I°(k,k), 
P 2h (k) = [ll(k)] 2 P^(k), 



(13) 
(14) 

(15) 



where the two terms represent contributions from two 
points in a single halo (lh) and points in different halos 
(2h) respectively. 

Likewise for the trispectrum, the contributions may be 
separated into those involving one to four halos 



T = T 



lh 



T 



2ft 



T 



3ft 



T 



lh 



(16) 



where here and below the argument of the trispectrum is 
understood to be (ki,k2, ka, kt). The term involving a 
single halo probes correlations of dark matter within that 
halo 

T lh = l2(ki, k 2 ,k 3 ,k 4: ) , (17) 

and is independent of configuration due to the assumed 
spherical symmetry for our halos. 

The term involving two halos can be further broken up 
into two parts 

T 2h = T 2h + T 2h^ ( 18 ) 
which represent taking three or two points in the first halo 



r 3 2 1 ' i = pp t (fc 1 )/ 3 1 (fc 2 ,A ; 3,fc4)/ 1 1 (fci; 
T 2 f = pp t (fc 12 )/ 2 1 (fc 1 , k 2 )4(h, h 



3 Perm., (19) 
- 2 Perm. (20) 
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Fig. 2. — The ratio of the single halo term contribution to that 
for a concentration width a\ n c — » for the (a) power spectrum and 
(b) trispectrum. The small scale behavior is increasingly sensitive 
to the high concentration tails for the higher order statistics. 



The permutations involve the 3 other choices of fej for the 
l\ term in the first equation and the two other pairings 
of the fcj's for the l\ terms in the second. Here, we have 
defined ki 2 = ki + k 2 ; note that k\ 2 is the length of one 
of the diagonals in the configuration. 

The term containing three halos can only arise with two 
points in one halo and one in each of the others 

T 3h = B pt( kl;k25 k 34 ) / i(fc3 j A :4 )/ 1 1 (fc 1 )A 1 (fc 2 ) (21) 
+pP t (k 1 )PP t (k 2 )ll(k 3} k 4 )ll(k 1 )ll(k 2 ) + 5 Perm. , 

where the permutations represent the unique pairings of 
the fti's in the I 2 factors. This term also depends on the 
configuration. The bispectrum in perturbation theory is 
given byf] 

^(k^k^k,.) = 2F 2 B (k p , k q )P(k p )P(k q ) + 2 Perm.^22) 



The kerne ls F s are Hpr ivpd in 



Goroff et al 



(1986) (see, equations 
A2 and A3 of Goroff et a] 1986; note that their P n = F n ), and we 
have written such that the symmetric form of -F n 's are used. TJhp 
use of the symmetric forrruaccounts for the factor of 2 in Eqs. (E2J) 
and factors of 4 and 6 in (|24|) . 
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Fig. 3. — Q sq at present broken into individual contributions un- 
der the halo description. The hierarchical model predicts a constant 
value for Q Bq in the deeply non-linear regime for clustering (HEPT). 
In the linear regime, the perturbation theory (PT) prediction is re- 
produced by the 4 halo term which is only ~ 1/2 of the total. See 
text for a discussion of discrepancies. 



with f| term given by second order gravitational pertur- 
bation calculations (see, below). 

Finally for four halos, the contribution is 



T 4h = j 



1 1 (fc 1 )/ 1 1 (fc2)/ 1 1 (fc 3 )/ 1 1 (fc4){TP t 

xP pt (k l )P pt (k 2 )P pt (k 3 ) + 3 Perm 



}, (23) 



where the permutations represent the choice of ki in the 
P's in the brac kets. The perturbation trispectrum can be 
written as (Fry 1984) 



T p = 

4[F|(k 12 ,-k 1 )F 2 s (k 12 ,k3)F(fc 1 )F(fc 12 )P(fc 3 ) + Perm.] 
+6[ J F 1 3 s (k 1 ,k 2 ,k3)F(fc 1 )P(fc 2 )F(fc 3 )+Pcrm.] . (24) 

The permutations involve a total of 12 terms in the first 
set and 4 terms in the second set. We now discuss the 
results from this modeling for a specific choice of halo input 
parameters and cosmology. 

2.3. Results 

2.3.1. Fiducial Model 

We evaluate the trispectrum under the halo model of 
the last sectio n assuming an NFW profile for the halos 
( Navarro et a] 1996) which depends on their virial mass M 
and concentration c. For the differential number density 
we take 



dn 
dMdc 

p{c)dc - 



( dn\ 

\dMJ ps 
1 



P(c) , 



(25) 



exp 



(In c — lnc) 

2<c 



dhic, 



where PS denotes the Prcss-Sch echter mass function. From 
the simulations of Bullock et a\ (2000), the mean and width 



of the concentration distribution is taken to be 

-0.13 



c(M, z) = 9(l + z)" 1 

fine = 0.2 , 



M 



M*(z) 



(26) 
(27) 



where M*(z) is the non- linear mass scale at which the 
peak-height threshold, f(M, z) = 1. 



Hi; 



This prescription differs from that in Cooray 
(2000) where cr\ nc — > since a finite distribution becomes 
increasingly important for the higher moments. To main- 
tain consistency we have also taken the mean concentra- 
tion directly from simulations rather than empirically ad- 
just it to match the power spectrum. For the same reason 
we choose a ACDM cosmological model with tt m = 0.3, 
Q\ = 0.7, h = 0.65 and a scale invariant spectrum of pri- 
mordial fluctuations. This model has mass fluctuations on 
the 8 h Mpc" 1 scale of ~ 1.0, consistent with the abun- 



dancc of galaxy clusters ( Viana & Liddle 1999) and COBE 
(Bunn & White 1997). For the linear power spectrum, we 
take the fitting f ormula for the transfer function given in 



Eiscnstein fc Hu| (1999). 



2.3.2. Comparisons 

In Fig. Ha), we show the logarithmic power spectrum 
A 2 (fc) = f?P(k)/2n 2 with contributions broken down to 
the lh and 2h terms today and the lh term at redshift of 
1. We find that there is an slight overprediction of power 
at scales corresp onding to 1 <; k <J 10 h Mpc" 1 when 
compared to the Peacock fc Dodds (1996) fitting function 
shown for redshifts of and 1, and a more substantial 
underprediction at small scales with ^ 10 h Mpc" 1 . 
Since the non-linear power spectrum has only been prop- 
erly studied out to overdensities A 2 ~ 10 3 with numerical 
simulations it is unclear whether the small-scale disagree- 
ment is significant. Fortunately it is on sufficiently small 
scales so as not to affect the lensing observables. 

For the trispectrum, we are mainly interested in terms 
involving T(ki, — ki, k 2 , — k 2 ), i.e. parallelograms which 
arc defined by either the length fci 2 or the angle between ki 
and k 2 . For illustration purposes we will take k\ = fc 2 and 
the angle to be 90° (k 2 — kj_) such that the parallelogram 
is a square. It is then convenient to define 



A s 2 q (fc) = 



2^2 



T 1 /3(k,-k,kj.,-k J 



(28) 



such that this quantity scales roughly as the logarithmic 
power spectrum itself A 2 (fc). This spectrum is shown in 
Fig. 0(b) with the individual contributions from the lh, 
2h, 3h, 4h terms shown. We test the sensitivity of our cal- 
culations to the width of the distribution in Fig. where 
we show the ratio between single halo contribution, as a 
function of the concentration distribution width, to the 
halo term with a delta function distribution a c = 0. As in 
the power spectrum the effect of increasing the width is to 
increase the amplitude at small scales due to the high con- 
centration tail of the distribution. Notice that the width 
effect is stronger in the trispectrum than the power spec- 
trum since the tails of the distribution are weighted more 
heavily in higher point statistics. 

To compare the specific scaling predicted by perturba- 
tion theory in the linear regime and the hierarchical ansatz 
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Table 1 



Dark Matter Power Spectrum Correlations 



k 


0.031 


0.044 


0.058 


0.074 


0.093 


0.110 


0.138 


0.169 


0.206 


0.254 


0.313 


0.385 


0.031 


1.000 


0.019 


0.041 


0.065 


0.086 


0.113 


0.149 


0.172 


0.186 


0.186 


0.172 


0.155 


0.044 


(-0.017) 


1.000 


0.036 


0.075 


0.111 


0.153 


0.204 


0.238 


0.261 


0.264 


0.251 


0.230 


0.058 


(0.023) 


(0.001) 


1.000 


0.062 


0.118 


0.183 


0.255 


0.302 


0.334 


0.341 


0.328 


0.305 


0.074 


(0.024) 


(0.024) 


(0.041) 


1.000 


0.102 


0.189 


0.299 


0.368 


0.412 


0.425 


0.412 


0.389 


0.093 


(0.042) 


(0.056) 


(0.027) 


(0.079) 


1.000 


0.160 


0.295 


0.404 


0.466 


0.485 


0.475 


0.453 


0.110 


(0.154) 


(0.076) 


(0.086) 


(0.094) 


(0.028) 


1.000 


0.277 


0.433 


0.541 


0.576 


0.570 


0.549 


0.138 


(0.176) 


(0.118) 


(0.149) 


(0.202) 


(0.085) 


(0.205) 


1.000 


0.434 


0.580 


0.693 


0.698 


0.680 


0.169 


(0.188) 


(0.180) 


(0.138) 


(0.229) 


(0.177) 


(0.251) 


(0.281) 


1.000 


0.592 


0.737 


0.778 


0.766 


0.206 


(0.224) 


(0.165) 


(0.177) 


(0.322) 


(0.193) 


(0.314) 


(0.396) 


(0.484) 


1.000 


0.748 


0.839 


0.848 


0.254 


(0.264) 


(0.228) 


(0.206) 


(0.343) 


(0.261) 


(0.355) 


(0.488) 


(0.606) 


(0.654) 


1.000 


0.858 


0.896 


0.313 


(0.265) 


(0.234) 


(0.202) 


(0.374) 


(0.259) 


(0.397) 


(0.506) 


(0.618) 


(0.720) 


(0.816) 


1.000 


0.914 


0.385 


(0.270) 


(0.227) 


(0.205) 


(0.391) 


(0.262) 


(0.374) 


(0.508) 


(0.633) 


(0.733) 


(0.835) 


(0.902) 


1.000 


\/% 


1.00 


1.01 


1.02 


1.03 


1.04 


1.07 


1.14 


1.23 


1.38 


1.61 


1.90 


2.26 



NOTES. — Diagonal normalized covariance matrix of the binned dark matter density field power spectrum with k values in units of h Mpc — . 
Upper triangle displays t he rrnm.ria.nrp found under the halo model. Lower triangle (parenthetical numbers) displays the covariance found in 
numerical simulations by [Meiksin fc White) (1999). Final line shows the fractional increase in the errors (root diagonal covariance) due to 
non-Gaussianity as calculated under the halo model. 



in the deeply non-linear regime, it is useful to define the 
quantity 



= T(k,-k, kL.-kj.) 

~ [8P 2 {k)P{V2k)}[4P 3 (k)} 



(29) 



In the halo prescription, Q sq at k ;> 10fc non iin ~ lOft-Mpc -1 
arises mainly from the single halo term. In perturbation 
theory Q sq w 0.085. The Q sq does not approach the per- 
turbation theory prediction as k — > since that contribu- 
tion appears only as one term in the 4 halo piece. There 
is an intrinsic shot noise error introduced by modeling the 
continuous density field by discrete objects. This error ap- 
pears large in the Q sq statistic since we have subtracted 
out the much larger connected (Gaussian) piece of the four 
point function. For example in the power spectrum covari- 
ance, the error induced by this approximation is much less 
than the Gaussian variance. 

The hierarchical ansatz predicts that Q sq = const, in 
the deeply non-linear regime. Its value is unspecified by 
the ansatz but is given as 



-jsat 
/ sq 



54 - 27 • 2 n + 2 • 3™ + 6" 
1 + 6 • 2 n + 3 • 3 n + 6 • 6" 



(30) 



under hyperexten ded perturbation theory (HEPT; pcocci 
marrc & Frieman). Here n — n(k) is the linear power spec- 
tral index at k. As shown in Fig. 0, the halo model predicts 
Q S q increases at high k. This behavior, also present at the 
three point level for the dark matter density field bispec- 
trum, suggests disagreement between the halo a pproach 
and hierarchical clustering ansatz (see, Ma & Fry 2000b), 
though numerical simulations do not yet have enough res- 
olution to test this disagreement. Fortunately the discrep- 
ancy is also outside of the regime important for lensing. 

To further test the accuracy of our halo trispectrum, we 
compare dark matter correlations predic ted by our method 
to those from numerical simulations by Meiksin fc White] 



(1999). For this purpose, we calculate the covariance ma- 
trix Cij from Eqn. (|J) with the bins centered at fcj and 
volume V a i = Airkfdki corresponding to their scheme. We 
also employ the parameters of their ACDM cosmology and 
assume that the parameters that defined the halo concen- 
tration properties from our fiducial ACDM model holds 
for this cosmological model also. The physical differences 
between the two cosmological model are minor, though 
normalization differences can lead to large changes in the 
correlation coefficients. 

In Table 0, we compare the predictions for the correla- 
tion coefficients 



Cij — 



C ; 



(31) 



33 



with the simulations. Agreement in the off diagonal el- 
ements is typically better than ±0.1, even in the region 
where non-Gaussian effects dominate, and the qualitative 
features such as the increase in correlations across the non- 
linear scale are preserved. 

A further test on the accuracy of the halo approach is to 
consider higher or der real-space m oments such as skewness 
and kurtosis. In Cooray & Hu (2000), we discussed the 
weak lensing convergence skewness under the halo model 
and f ound it to be in agreement with numerical predictions 
from IWhite fc Hu| (1999). The fourth moment of the den- 



sity field, under certain approxim ations, was calculated by 
Scoccimarro, Zaldarriaga & Hui (1999) using dark matter 
halos and was found to be in good agreement with N- 
body simulations. Given that density field moments have 



already been studied by Scoccimarro, Zaldarriaga & Hui 
we no longer consider them here other than to suggest 
that the halo model has provided, at least qualitatively, a 
consistent description better than any of the perturbation 
theory arguments. 

Even though the dark matter halo formalism provides a 
physically motivated means of calculating the statistics of 
the dark matter density field, and especially higher order 




Fig. 4. — Weak lensing convergence (a) power spectrum and (b) 
trispectrum under the halo description. Also shown in (a) is the 
prediction from the PD nonlinear power spectrum fitting function. 
We have separated individual contributions under the halo approach 
and have assumed that all sources are at z B = 1 . We have also shown 
the shot noise contribution to the power spectrum assuming a survey 
down to a limiting magnitude of R ~ 25 with an intrinsic rms shear 
of 0.4 in each component. 



correlations, there are several limitations of the approach 
that should be borne in mind when interpreting results. 
The approach assumes all halos to share a parameterized 
spherically-symmetric profile. We have attempted to in- 
clude variations in the halo profiles with the addition of 
a distribution function for concentration parameter based 
on results from numerical simulations. Unli ke our pre- 



vious calculat ions presented in Cooray et a] (2000) and 
Cooray & Hrj (2000), we have not modified concentration- 
mass relation to fit the PD non-linear power spectrum, 
but rather have taken results directly from simulations as 
inputs. Though we have partly accounted for halo profile 
variations, the assumption that halos are spherical is likely 
to affect detailed results on the configuration dependence 
of the trispectrum. Since we are considering a weighted 
average of configurations, our tests here are insufficient to 
establish the validity of the trispectrum modeling in gen- 
eral. Further numerical work is required to quantify to 



what extent the present approach reproduces simulation 
results for the full trispectrum. 

3. CONVERGENCE POWER SPECTRUM COVARIANCE 
3.1. General Definitions 

Weak lensing probes the statistical properties of the 
shear field on the sky which is a weighted projection of 
the matter distribution along the line of sight to the source 
galaxies. As such, the observables may be reexpressed as 
a scalar quantity, the convergence k, on the sky. 

Its power spectrum and trispectrum are defined in the 
flat sky approximation in the usual way 

<«(l 1 ) K (l 2 )) = (27r) 2 <5 D (l 12 )Cr, 
(k(1i) ■ • • k(U)) c = (27r) 2 fc(l 1234 )T K (l 1 , 1 2 , 1 3 , 1 4 ) . (32) 

These are related to the den sity po wer s pectrum and 



trispectrum by the projections (Kaiser 1992; Bcoccimarro, 
Zaldarriaga fc Hui| 1999) 



rpK, 



W(r) 2 

d A 

dr^T 

4 



i 

d A ' 

h_ J2_ hi_ I4 



(33) 
(34) 



where r is the comoving distance and d A is the angular 
diameter distance. When all background sources are at a 



10" 



< 



10' 



(a) 




Gaussian 

with nonGaussianities (total) 
with nonGaussianities (1 h term ) 



RlO" 1 

o 

< 



10" 



10' 




Fig. 5. — The fractional errors in the measurements of the con- 
vergence band powers. In (a), we show the fractional errors under 
the Gaussian approximation, the full halo description, the Gaussian 
plus single halo term, and the Gaussian plus shot noise term (see 
y.3[ ). As shown, the additional variance can be modeled with the 
single halo piece while shot noise generally becomes dominant be- 
fore non-Gaussian effects hprnmp largp Tn (U\ we compare the halo 
model with simulations from White & Hu (1999). The decrease in 
the variance at small I in the simulations is due to the conversion of 
variance to covariance by the finite box size of the simulations. 
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Table 2 

Weak Lensing Convergence Power Spectrum Correlations 





97 


138 


194 


271 


378 


529 


739 


1031 


1440 


2012 


97 


1.00 


0.04 


0.05 


0.07 


0.08 


0.09 


0.09 


0.09 


0.08 


0.08 


138 


(0.26) 


1.00 


0.08 


0.10 


0.11 


0.12 


0.12 


0.12 


0.11 


0.11 


194 


(0.12) 


(0.31) 


1.00 


0.14 


0.17 


0.18 


0.18 


0.17 


0.16 


0.15 


271 


(0.10) 


(0.21) 


(0.26) 


1.00 


0.24 


0.25 


0.25 


0.24 


0.22 


0.21 


378 


(0.02) 


(0.09) 


(0.24) 


(0.38) 


1.00 


0.33 


0.33 


0.32 


0.30 


0.28 


529 


(0.10) 


(0.14) 


(0.28) 


(0.33) 


(0.45) 


1.00 


0.42 


0.40 


0.37 


0.35 


739 


(0.12) 


(0.16) 


(0.17) 


(0.34) 


(0.38) 


(0.50) 


1.00 


0.48 


0.45 


0.42 


1031 


(0.15) 


(0.18) 


(0.15) 


(0.27) 


(0.33) 


(0.48) 


(0.54) 


1.00 


0.52 


0.48 


1440 


(0.18) 


(0.15) 


(0.19) 


(0.19) 


(0.32) 


(0.36) 


(0.53) 


(0.57) 


1.00 


0.54 


2012 


(0.19) 


(0.22) 


(0.16) 


(0.32) 


(0.27) 


(0.46) 


(0.50) 


(0.61) 


(0.65) 


1.00 



NOTES. — Covariance of the binned power spectrum when sources are at a redshift of 1. Upper triangle displays 
the halo model. Lower triangle (parenthetical numbers) displays the covariance found in numerical simulations by White 
consistent with these simulations, we use the same binning scheme as the one used there. 



hp cnvarianr e found under 



Hu (1999). To be 



distance of r s , the weight function becomes 



Wir) 



H 2 d A {r)d A (r s -r) 
c 2 a d A (r s ) 



(35) 



for simplicity, we will assume r s = r(z s = 1). In deriving 
Eq. (p4|), we have used the Limber approximation (Limbci 
1954)by setting k — l/d A and the flat-sky approximation. 
A potential problem in using the Limber approximation is 
that we implicitly integrate over the unperturbed photon 
paths (Born approximation) . The Born a pproximat ion has 
been tested in numerical simulations by Jain ct a| (2000; 
see their Fig. 7) and found to be an excellent approxima- 
tion for the two point statistics. The same approxima- 
tion can also be tested through lens-lens coupling involv- 
ing lenses at two different redshifts. For higher order cor- 
relations, analytical calculation s in the mildl y non-linear 
regime by Van Waerbeke et &\ (2000b; also, Bernardeau 
et alQ997; [Schneider ct al| 1998) indicate that corrections 
are again less than a few percent. Thus, our use of the 
Limber approximation by ignoring the lens-lens coupling 
is not expected to change the final results significantly. 

For the purpose of this calculation, we assume that 
upcoming weak lensing convergence power spectrum will 
measure binned logarithmic band powers at several It's in 
multipole space with bins of thickness 5k. 



f d 2 l I 2 , 
J si A si 2w 



(36) 



where A s (li) — f d 2 l is the area of 2D shell in multipole 
and can be written as A a (k) = 2wli5li + Tr(Sli) 2 . 
We can now write the signal covariance matrix as 



Ci 



rpK 

ij 



1 

A 



(2n) 2 2 



A s i 



d 2 h 



l 2 l 2 



A sj (2tt)< 



T K (l- -1 1 -1 



(37) 
(38) 



where A is the area of the survey in steradians. Again 
the first term is the Gaussian contribution to the sample 
variance and the second the non-Gaussian contribution. A 



realistic survey will also have shot noise variance due to 
the finite number of source gala xies in the survey. We will 



return to this point in the §3.3 



3.2. Comparisons 

Using the halo model, we can now calculate contribu- 
tions to lensing convergence power spectrum and trispec- 
trum. The logarithmic power spectrum, shown in Fig. 0(a) , 
shows the same behavior as the density field when com- 
pared with the PD results: a slight overprediction of power 
when I > 10 3 . However, these differences are not likely to 
be observable given the shot noise from the finite number 
of galaxies at small scales. 

In Fig |4|(b) , we show the scaled trispectrum 



A« q (0 = ^T"(l,-l,l x ,-lx) 1/3 . 

27T 



Table 3 

Inverse Fisher Matrix (x 10 3 ) 



(39) 



Pi 


A 


In A 


n K 


n s 




n A 


1.57 


-5.96 


-1.39 


4.41 


-1.76 


In A 




25.89 


5.83 


-17.34 


6.74 


n K 






1.41 


-3.81 


1.43 


n s 








14.01 


-6.03 












2.67 




Pi 


A 


In A 


n x 


n s 




n A 


2.03 


-7.84 


-1.82 


5.76 


-2.30 


In A 




33.92 


7.65 


-22.79 


8.91 


n K 






1.78 


-5.01 


1.95 


n s 








18.43 


-7.85 












3.44 



NOTES. — Inverse Fisher matrix under the Gaussian assumption 
(top) and the halo model (bottom). The error on an individual 
parameter is the square root of the diagonal element of the Fisher 
matrix for the parameter while off-diagonal entries of the inverse 
Fisher matrix shows correlations, and, thus, degeneracies, between 
parameters. We have assumed, a full sky survey (/ s k y = 1) with 
parameters as described in 



8 



10" 



o 
O 



1 
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: 






2802 






5432 






10540 






20464 






39710 




/ 1h 


77072 


/ 1 






Fig. 6. — (a) The correlation coefficient, 
multipolc h with lj as shown in the figure. We show the correlations 
calculated with the full halo model and also with only the single halo 
term for lj = 77072. In (b), we show the non-Gaussian correlation 
coefficient Cf^^ , which only involves the trispectrum (see, Eq. ^). 
The transition to full correlation is due to the domination of the 
single halo contribution. 



White Sz Hu (1999) simulations. The decrease in errors 
from the simulations at small I reflects finite box effects 
that convert variance to covariance as the fundamental 
mode in the box becomes comparable to the bandwidth. 
The correlation between the bands is given by 



Cij = 



Cs. 



(41) 



j j 



In Table B we compare the halo predictions to the simula- 
tions by |White fc Hu (1999). The upper triangle here is 



the correlations under the halo approach, while the lower 
triangle shows the correlations found in numerical simula- 
tions. The correlations along individual columns increase 
(as one goes to large l's or small angular scales) consistent 
with simulations. In Fig. we show the correlation coeffi- 
cients with (a) and without (b) the Gaussian contribution 
to the diagonal. 

We show in Fig. ^(a) the behavior of the correlation co- 
efficient between a fixed lj as a function of U. When U = lj 
the coefficient is 1 by definition. Due to the presence of 
the dominant Gaussian contribution at U = lj, the coef- 



L and 



j_i that decreases as lj increases and non-Gaussian 



ficient has an apparent discontinuity between l l — tJ 
k =1 

effects dominate 

To better understand this behavior it is useful to isolate 
the purely non-Gaussian correlation coefficient 



C, 



NG 



Ti, 



yjTuTij 



(42) 



As shown in Fig. ^](b), the coefficient remains constant for 
hi <C lj and smoothly increases to unity across a transition 
scale that is related to where the single halo terms starts to 
contribute. A comparison ofFig.|(b) and|(b), shows that 
this transition happens around / of few hundred to 1000. 
Once the power spectrum is dominated by correlations in 
single halos, the fixed profile of the halos will correlate 
the power in all the modes. The multiple halo terms on 
the other hand correlate linear and non-linear scales but 
at a level that is generally negligible compared with the 
Gaussian variance. 

The behavior seen in the halo based covariance, how- 
ever, is not present when the covariance is calculat ed with 



wher q lj_ — I and 1 ■ lj^ — 0. The projected lenaing tri3pcc - 
trum l again shows the same behavior as the density field 



hierarchical arguments fo r the trispectrum (see, Scocci- 



marro, Zaldarriaga fc Huj 1999). With hierarchical argu- 



trispectrum with similar conditions on k^ 's. 

We can now use this trispectrum to study the contri- 
butions to the covariance, which is what we are primarily 
concerned here. In Fig. |^a, we show the fractional error, 



Ad 

Ci 



C, 



(40) 



for ban ds lj given in T able following the binning scheme 
used by |White fc Hu[ (1999) on 6° x 6° fields. The dashed 
line compares that with the Gaussian errors, involving the 
first term in the covariance (Eq. [58]). At multipoles of a 
few hundred and greater, the non-Gaussian term begins 
to dominate the contributions. For this reason, the errors 
are well approximated by simply taking the Gaussian and 
single halo contributions. 

In Fig. ||b, we compare these results with those of the 



ments, which are by construction only valid in the deeply 
nonlinear regime, one predicts correlations which are, in 
general, constant across all scales and shows no decrease 
in correlations between very small and very large scales. 
Such hierarchical models also violate the Schwarz inequal- 
ity with correlati ons greater than 1 between larg e and 
small scal es (e.g., Scoccimarro, Zaldarriaga & Hui 1999; 
Hamilton 2000). The halo model, however, shows a de- 
crease in correlations similar to numerical simulations sug- 
gesting that the halo model, at least qualitatively, provides 
a better approach to studing non-Gaussian correlations in 
the translinear regime. 

3.3. Effect on Parameter Estimation 

Modeling or measuring the covariance matrix of the 
power spectrum estimates will be essential for interpret- 
ing observational results. In the absence of many fields 
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where the covariance can be estimated directly from the 
data, the halo model provides a useful, albeit model de- 
pendent, quantification of the covariance. As a practi- 
cal approach one could imagine taking the variances esti- 
mated from the survey under a Gaussian approximation, 
but which ac counts for uneven sampling and edge effects 
( Hu fc White 2000), and scaling it up by the non-Gaussian 



to Gaussian variance ratio of the halo model along with in- 
clusion of the band power correlations. Additionally, it is 
in principle possible to use the expected correlations from 
the h alo model to decorrelate individual band power mea- 



surem ents, similar to studies involving CMB temperature 



anisotropy and galaxy power spectra (e.g., Hamilton 1997; 
Hamilton k Tegmark| 2000). 



We can estimate the resulting effects on cosmological 
parameter estimat ion with an ana logous procedure on the 
Fisher matrix. In Hu & Tegmark (1999), the potential of 
wide-field lensing surveys to measure cosmological parame- 
ters was investigated using the Gaussian approximation of 
a diagonal covariance and Fisher matrix techniques. The 
Fisher matrix is simply a projection of the covariance ma- 
trix onto the basis of cosmological parameters pi 



Fq/3 = 



dp/3 ' 



(43) 



where the total covariance includes both the signal and 
noise covariance. Under the approximation of Gaussian 
shot noise, this reduces to replacing C ; K — > C* + Cf N in 
the expressions leading up to the covariance equation (p8|), 



The stint, nnisp pnwpr spppt.mm is givpn hy 



qSN _ (lint) 



(44) 



where (71, 



a/2 



0.4 is the rms noise per component in- 



troduced by intrinsic ellipticities and measurement errors 
and n ~ 6.6 x 10 8 sr _1 is the surface number density of 
background source galaxies. The numerical values here are 
appropriate f or surveys th at reach a limiting magnitude in 
R ~ 25 (e.g., |5mail ct a\ 1995). 

Under the approximation that there are a sufficient 
number of modes in the band powers that the distribution 
of power spectrum estimates is approximately Gaussian, 
the Fisher matrix quantifies the best possible errors on 
cosmological parameters that can be achieved by a given 
survey. In particular F~ x is the optimal covariance matrix 

of the parameters and (-F -1 )*/ 2 is the optimal error on the 
ith parameter. Implicit in this approximation of the Fisher 
matrix is the neglect of information from the cosmologi- 
cal parameter dependence of the covariance matrix of the 
band powers themselves. Since the covariance is much less 
than the mean power, we expect this information content 
to be small. 

In order to estimate the effect of non-Gaussianities on 
the cosmological parameters, we calculate the Fisher ma- 
trix elements using our fiducial ACDM cosmological model 
and define the dark matter density field, today, as 



A 2 (fc) = ^( — 



n a +3 



T 2 (k) . 



(45) 



Here, A is the amplitude of the present day density fluc- 
tuations and n s is the tilt at the Hubble scale. The den- 
sity power spectrum is evolved to higher redshifts using 



the growth function G(z) ( Peebles 1980) and the transfer 
function T(k) is calculated using the fitting functions from 



Eiscnstein & Hu (1999). Since we are only interested in the 



relative effect of non-Gaussianities, we restrict ourselves to 
a s mall subset of th e cosmological parameters considered 
by Hu fc Tegmarkj (1999) and assume a full sky survey 



with / sky = 1 

In Table 3, we show the inverse Fisher matrices deter- 
mined under the Gaussian and non-Gaussian covariances, 
respectively. For the purpose of this calculation, we adopt 
the bi nning scheme as shown in Table 2, following White 
fc Hu| (1999). The Gaussian errors are computed using the 
same scheme by setting T K — 0. As shown in Table 3, the 
inclusion of non-Gaussianities lead to an increase in the 
inverse Fisher matrix elements. We compare the errors 

on individual parameters, mainly {F^ 1 ) 1 ^ 2 , between the 
Gaussian and non-Gaussian assumptions in Table ^[ The 
errors increase typically by ~ 15%. Note also that band 
power correlations do not necessarily increase cosmological 
parameter errors. Correlations induced by non- linear grav- 
ity introduce larger errors in the overall amplitude of the 
power spectrum measuremenents but have a much smaller 
effect on those parameters controlling the shape of the 
spectrum. 

For a survey of this assumed depth, the shot noise power 
becomes the dominant error before the non-Gaussian sig- 
nal effects dominate over the Gaussian ones. For a deeper 
survey with better imaging, such as the one planne d with 



Large-ap erture Synoptic Survey Telescope (LSST; Tyson 
fc Angel| 2000)0, the effect of shot noise decreases and non- 



Gaussianity is potentially more important. However, the 
non-Gaussianity itself also decreases with survey depth, 
and as we now discuss, in terms of the effect of non- 
Gaussianities, deeper surveys should be preferred over the 
shallow ones. 

3.4. Scaling Relations 

To better understand how the non-Gaussian contribu- 
tion scale with our assumptions, we consider the ratio of 
i to the Gi 
Hul| 1999), 



non-Gaussian variance to the Gaussian variance (Scocci- 
marro, Zaldarriaga 



with 



_ A si T£ 



(2tt) 2 2C| 



(46) 



(47) 



Under the assumption that contributions to lensing con- 
vergence can be written through an effective distance r*, at 
half the angular diameter distance to background sources, 
and a width Ar for the lensing window function, the ratio 
of lensing convergence trispectrum and power spectrum 
contribution to the variance can be further simplified to 



R 



U off P 2 (r*) 



(48) 



Since the lensing window function peaks at r*, we have 
replaced the integral over the window function of the den- 
sity field trispectrum and power spectrum by its value at 

3 http:/ /www. dmtelescope.org 
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the peak. This ratio shows how the relative contribution 
from non-Gaussianities scale with survey parameters: (a) 
increasing the bin size, through A S i (oc SI), leads to an 
increase in the non-Gaussian contribution linearly, (b) in- 
creasing the source redshift, through the effective volume 
of lenses in the survey (V c g ~ r^Ar), decreases the non- 
Gaussian contribution, while (c) the growth of the density 
field trispectrum and power spectrum, through the ratio 
T/P 2 , decreases the contribution as one moves to a higher 
redshift. The volume factor quantifies the number of fore- 
ground halos in the survey that effectively act as gravi- 
tational lenses for background sources; as the number of 
such halos is increased, the non-Gaussianities are reduced 
by the central limit theorem. 

In order to determine whether its the increase in volume 
or the decrease in the growth of structures that lead to a 
decrease in the relative importance of non-Gaussianities 
as one moves to a higher source redshift, we numerically 
calculated the non-Gaussian to Gaussian variance ratio un- 
der the halo model for several source redshifts and survey 
volumes. Up to source redshifts ~ 1.5, the increase in vol- 
ume decreases the non-Gaussian contribution significantly. 
When surveys are sensitive to sources at redshifts beyond 
1.5, the increase in volume becomes less significant and the 
decrease in the growth of structures begin to be important 
in decreasing the non-Gaussian contribution. Since, in the 
deeply non-linear regime, T/P 2 scales with redshift as the 
cube of the growth factor, this behavior is consistent with 
the overall redshift scaling of the volume and growth. 

Given that scalings always lead to decrease the effect of 
non-Gaussianities in deep lensing surveys, with a decrease 
in the shot noise contribution also, shallow surveys that 
only probe sources out to redshifts of a few tenths are 
more likely to be dominated by non-Gaussianities; such 
shallow surveys are also likely to be affected by i ntrin- 
sic correlations o f gala xy shapes (e.g. , |Catelan et a] 2000; 



Croft fc Metzler 2000; Heavens et al. 2000). These possi 



bilities, generally, suggests that deeper surveys should be 
preferred over shallow ones for weak lensing purposes. 

4. CONCLUSIONS 
Weak gravitational lensing due to large scale struc- 



ture rtroviclcs important mtormation on trie evolution ot 
clusteri ng and angular : 



rngmrllngiral parampt.prs TViis infnrmat.inn rnmplomonts 

what can be learned from cosmic microwave background 
anisotropy observations. The tremendous progress on the 
observational front warrants detailed studies of the statis- 



Table 4 
Parameter Errors 



In A 



-a* 



Gausslanl 



n n3Q n im n 037 n ns n nsi 



Full 

Increase (%) 



0.045 
15.3 



0.184 
15 



0.042 
13.5 



0.135 
14.4 



0.058 
13.7 



NOTES. — Parameter errors, (F 1 )^/ 2 , under the Gaussian assump- 
tion (top) and the halo model (bottom) and following the inverse- 
Fisher matrices in Table 3. We have assumed a full sky survey 
(/sky = 1) with parameters as described in § 3.3 



tical properties of the lensing observables and their use in 
constraining cosmological models. 

The non-linear growth of large-scale structure induces 
high order correlations in the derived shear and conver- 
gence fields. In this work, we have studied the four point 
correlations in the fields. Four point statistics are special 
in that they quantify the errors in the determination of 
the two point statistics. To interpret future lensing mea- 
surements on the power spectrum, it will be essential to 
have an accurate assessment of the correlation between the 
measurements . 

Using the halo model for clustering, we have provided 
a semi-analytical method to calculate the four point func- 
tion of the lensing convergence as well as the dark matter 
density field. We have tested this model against numeri- 
cal A^-body simulations of the power spectrum covariance 
in both the density and convergence fields and obtained 
good agreement. As such, this method provides a practical 
means of estimating the error matrix from future surveys 
in the absence of sufficiently large fields where it may be 
estimated directly from the data or large suites of A'-body 
simulations where it can be quantified in a given model 
context. Eventually a test of whether the covariance ma- 
trix estimated from the data and the theory agree may 
even provide further cosmological constraints. 

This method may also be used to study other aspects 
of the four point function in lensing and any field whose 
relation to the dark matter density field can be modeled. 
Given the approximate nature of these approximations, 
each potential use must be tested against simulations. 
Nonetheless, the halo model provides the most intuitive 
and extensible means to study non-Gaussianity in the cos- 
mological context currently known. 

We acknowledge useful discussions with Dragan Huterer, 
Roman Scoccimaro, Uros Scljak, Ravi Sheth and Matias 
Zaldarriaga. WH is supported by the Keck Foundation. 
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